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ABSTRACT 


The  Ildar  system  oiay  now  be  considered  as  being  essen- 
tially completed.  A detailed  description  in  its  present, 
final  state,  Including  data  acquisition  and  processing, 
is  given.  Most  recent  extensions  of  the  system  are  photon 
counting  with  range  gating  facilities,  as  well  as  electro- 
mechanical fluorescence  suppression,  which  both  provide 
a considerable  increase  in  sensitivity  for  long-range, 
stratospheric  aerosol  sensing  up  to  more  than  30  km  alti- 
tude. Examples  of  all  inodes  of  data  recording  including 
photon  counting  are  presented.  Furthermore,  examples  of 
simultaneous  lidar  backscatter  measurements  and  cable- 
car  recordings  of  profiles  of  Aitken  nuclei  concentration 
and  aerological  parameters  are  discussed,  A first  high- 
altitude  measurement  using  the  photon  counting  method  is 
presented. 

Extensive  mathematical  and  theoretical  procedures  have 
been  worked  out  for  calibration  of  the  system,  using 
simultaneously  acquired  aerosol  and  aerological  data,  and 
for  evaluation  of  the  return  signals  once  the  system  is 
absolutely  calibrated.  The  procedures  are  described  in 
detail,  examples  of  calibration  measurements  are  presented 
and  discussed. 
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OBJECTIVES  AND  BACKGROUND 


. Preliminary  Remarks 


The  present  report  subsequently  continues  and  extends 
information  given  in  the  preceding  Annual  Report  of  June 
1974  (1)  under  the  same  grant  (DA-ER0-591-73-GO057),  In 
order  to  avoid  repetition,  reference  with  regard  tc  the 
General  Background  is  made  to  (1). 


ytain  Objective  of  Work  in  Reporting  Period 


The  objective  of  our  work  under  above  grant  is  to  develop 
a calibrated  lidar  system  for  aerosol  x'emote  sensing,  and 
to  employ  it  for  studying  aerosol  transport,  especially 
for  monitoring  the  vertical  aerosol  distribution  as  func- 
tion of  serological  parameters. 

The  course  our  efforts  took  to  develop  the  lidar  system 
as  such,  as  well  as  to  establish  the  necessary  theoretical 
basis  for  utilizing  data  obtained,  has  revealed  that  con- 
siderable work  will  yet  have  to  be  expended  to  attain  the 
objective  set.  We  therefore  considered  it  expedient  during 
the  reporting  period  to  concentrate  our  resources  on  large- 
ly completing  two  branches  of  our  basic  development  work: 


a)  technological  development  of  the  lidar  system  and 
of  the  requisite  electronic  data  processing  system; 

b)  establishing  the  theoretical  basis  and  computer 
programs  for  numeric  evaluation  of  the  recorded 
data. 

By  contrast,  the  performance  of  series  of  actual  measure- 
ments has  been  temporarily  postponed.  We  have  limited  our 
work  to  collecting  lidar  data,  aerosol  parameters  obtained 
at  our  stations  and  by  the  cable-car  telemetry  system, 
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and  aerological  data  on  single  days  of  especially  favor- 
able atmospberic  conditions,  to  make  available  some 
material  for  the  initial  testing  of  the  developed  lidar 
calibration  method. 

Detailed  reports  on  the  above  items  (1)  and  (2)  are  found 
in  the  chapters  below.  The  functioning  of  the  lidar  system 
at  its  final  stage  of  completion,  and  of  the  evaluating 
system,  is  represented  in  some  examples.  Likewise  the 
lidar  calibration  method  developed  by  us,  is  presented  by 
some  numerical  examples. 

Completion  of  lidar  calibration  as  well  as  practical  employ- 
ment of  the  entire  system  for  aerosol  remote  sensing  are 
intended  to  be  accomplished  in  the  scope  of  the  next  grant. 

II.  OBTAINING  CF  DATA  JULY  1974  - JUNE  1973 

Those  facilities  of  our  Institute  that  are  essential  with- 
in the  scope  of  the  grant,  were  employed  partly  to  run  con- 
tinuous measurements,  partly  to  conduct  series  of  meas- 
urements or  single  measurements,  the  extent  of  which  is 
briefly  summarized  as  follows. 

1 . Rawinsonde 


A total  of  22  ascents  with  an  average  summit  height  of 
30  km,  were  conducted  and  evaluated,  partly  for  purposes 
of  lidar  calibration,  partly  within  the  scope  of  basic 
atmospheric  investigations. 

2.  Cable  Car  Telemetry  Data 

A total  of  1441  measuring  runs  were  conducted  by  the  Zug- 
spitze  cable-car  telemetry  system  in  the  reporting  period. 


and  were  largely  evaluated.  The  aerological  data,  and  pro 
files  of  aerosol  concentration  (Aitken  nuclei)  were  used, 
in  part,  for  comparison  against  the  lidar  backscatterlng 
profiles. 

3.  Measurements  Taken  at  Three  Stations 

Following  measurements  and  recordings  have  been  continued 
without  any  interruptions: 

a)  Natural  radioactivity  in  the  air  (RaB,  RaC); 

b)  Condensation  nuclei  with  size  distribution; 

c)  Impactor  measurements; 

d)  Atmospheric  electricity; 

e)  Meteorology. 

A.  Lidar  Backscatterlng  Profiles 

The  extent  of  the  lidar  backscatterlng  profiles  taken  is 
as  follows? 

For  purposes  of  lidar  calibration,  synchronous  measure- 
ments were  conducted  by  lidar,  rawinsonde  and  impactors, 
on  12  days. 


III.  FINAL  TECHillCAL  STATUS  OF  THE  LIDAR  SYSTEM,  INCLU- 

DING  DATA  PROCESSING 


The  lidar  system,  which  had  been  put  into  operation  in 
1973,  had  been  described  in  its  basic  version  in  the  5th 
annual  report.  During  the  following  two  years  the  system 
had  been  extended  by  a number  of  supplementary  facilities 
Now  it  may  be  considered  as  completed  and  shall  be  de- 
scribed as  a whole  in  detail.  The  mechanical  layout  of 
the  system  is  shown  in  Fig. ill. 1,  the  block  scheme 
Fig. III. 2 shows  the  electronics  in  detail. 


Laser  Transmitter  with  Power  Supply  (cf.  Figs  III.1 
and  III. 2) 


The  Q-switched  ruby  laser  transmitter  (1)  is  mounted  hori- 
2onta‘‘ly  cn  a vibration-free  frame  (2)  within  a penthouse 
on  the  roof  of  the  Institute  building. 

The  holographic  quality,  Czochralski  type  ruby  (length 
165  mm,  diameter  9.5  mm)  is  installed  in  a solid  metal 
block  (5)  with  an  elliptical,  silver  coated  cavity,  to- 
gether with  the  flash  lamp,  and  is  water-cooled.  The  res- 
onator of  approximately  1.500  mm  length  consists  of  a 
highly  reflecting  end  mirror  in  the  rear  (4)  and  a front 
output  mirror  (5)  (BK7  glass  plate,  reflection). 
Filtered  and  dried  air  from  an  Institute-own  compressed- 
air  system  is  blown  on  the  ruby  end  faces  to  protect 
them,  particularly  from  condensed  water.  Q-switching  is 
done  ty  a Kerr  cell  (8)  and  a polarizer  (7). 

The  output  power  is  100  Megawatts  meximum,  the  pulse 
length  about  20  nanoseconds.  This  results  in  a pulsf" 
energy  of  2 Joules.  The  beam  divergence  is  5 mrad  total 
(no  expanding  telescope  is  used). 

Behind  the  end  mirror,  a battery-operated  trigger  sensor 
(8)  is  mounted,  which  picks  up  the  residual  laser  light 
penetrating  the  mirror  by  means  of  a photodiode,  and  gen- 
erates a trigger  pulse  for  the  data  acquisition  elec- 
tronics . 

The  power  supply  for  the  flashlamp,  (9),  which  includes 
the  Kerr  cell  driving  unit  (10)  are  located  in  the  same 
room.  It  is  amply  designed  for  permitting  a pulse  repe- 
tition  frequency  cf  up  to  2 sec  at  the  full  output 
energy  of  4 kilowattseconds  per  pulse.  This  is  of  great 
benefit  in  observing  rapidly  changing  processes.  For  sta- 
bility reasons,  howevor,  the  repetition  frequency  is  to 
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be  limited  to  1 set “ * . Through  an  aperture  in  the  wall 
of  the  loom,  the  laser  bean  hits  a 45  deg.  mirror  and  is 
directed  upwards  (see  below). 

2.  Frequency  Doubler 


A KD*P  crystal  (11)  is  used  for  geneiating  the  doubled 
ruby  Crequency,  347  nanometers.  A filter  cuvette  (12)  con- 
taining copper  sulphate  solution  removes  the  residual  un- 
doubled laser  light.  Both  crystal  and  cuvette  are  mounted, 
cogether  with  one  45  degree  dielectric  first  surface 
mirrcr  (13)  for  each  transmitter  frequency,  on  a carriage 
(14)  just  above  the  receiving  telescope  secondary  mirror 
(descx’ibed  later).  By  means  of  a servo  motor  operated  from 
the  operating  console  (15)  the  carriage  con  be  moved  from 
one  position  to  the  other,  thus  swinging  the  doubling 
equipment  ^nto  and  out  of  the  laser  beam  path  and  properly 
positioning  either  deviating  mirror.  Since  the  doubling 
efficiency  of  the  KD*P  crystal  (about  8%  maximum)  depends 
very  critically  on  the  exact  positioning  of  its  optical 
sucis  with  respect  to  the  laser  beam  (index  matching),  it 
is  alignabie  0/  two  precision  screws.  The  switching  time 
is  about  20  seconds,  allowing  double  frequency  firings 
within  short  time  intervals. 

3.  Receiver  Unit 


3.1.  Telescope 

The  backscattered  laser  light  is  picked  up  by  an  astro- 
nomical quality,  52  cm  dia.  Cassegrain  reflecting  tele- 
scope (16),  consisting  of  a concave  primary  mirror  (17) 
with  a borehole  in  its  center  and  a convex  secondary 
mirror  (IS)  generating  a secondary  focus  just  behind  the 
borehole.  The  resulting  equivalent  focal  length  is  3.8  m. 
The  telescope  is  mounted  in  a fixed  vertical  position 
within  a barrel-shaped  waterproof  case  (19)  permanently 
installed  in  the  open  on  the  observation  platform.  It  is 


closed  on  top  by  a tilted  precision  quartz  glass  plate 
allowing  rain  and  snow  water  to  flow  off  easily,  with  the 
transmitter  exit  window  (20)  placed  horizontally  in  its 
middle.  This  window  also  consists  of  quartz  glass  provi- 
ding good  transmission  for  the  3^7  nm  wavelength. 

The  laser  beam  enters  horizontally  through  a lateral 
flemge  and  hits  the  deviating  mirror/frequency  doubler 
assembly  mentioned  above,  and  leaves  the  case  through  the 
exit  window.  Despite  some  disadvantages,  e.g.  stray  light 
generation  in  the  case,  we  chose  this  concentric  con- 
figuration to  be  able  to  receive  returns  from  all  rauiges 
of  interest  without  any  readjusting. 

A temperature  regulated  heating  prevents  water  condensate 
formation  and  keeps  the  temperature  of  the  frequency 
doubling  crystal  constant.  This  is  necessary  because  the 
optimal  position  of  the  crystal  depends  on  temperature. 

Within  the  focal  plane  of  the  telescope  an  iris  stop  (21) 
limiting  the  field  of  view  of  the  system,  and  an  optical 
shutter  (22)  operated  electrically  from  the  console  are 
mounted.  Below  the  shutter,  a filter  turret  (23)  containing 
a 694  and  a 347  nm  interference  filter,  bandwith  3 nm, 
for  background  reduction,  and  a gray  wedge  (24)  are  placed. 
Both  units  are  also  operated  from  the  console  by  magslips. 
By  means  of  the  gray  wedge,  the  sensitivity  of  the  re- 
ceiver can  be  reduced  if  necessary. 

3.2.  Photomultiplier  Tube  Assembly 


The  signals  are  received  by  an  EMI  9816  photomultiplier 
tube  (PMT)  with  S 20  cathode  (25),  which  is  located  in- 
side a commercial  cooling  box  (26).. In  order  to  reduce 
the  dark  current  the  PMT  can  be  cooled  down  to  minus 
20  deg.  by  means  of  Peltier  elements.  The  light  enters 
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3.3.  t -Control 


the  cooling  box  through  a double-walled,  heated  quartz 
window  (27)  which  is  evacuated  to  prevent  condensate 
formation.  Outside  the  cooling  box  u collimating  lens 
(28)  and  a hjlder  for  a polarizer  (29)  is  installed.  By 
means  of  the  polarizer  the  degree  of  polarization  of  the 
backscattered  light  can  be  measured.  For  this  purpose  the 
box  is  mounted  on  a rotatable  support  with  a vernier 
scale  (30).  With  the  onset  of  double  frequency  operation 
polarization  measurements  have  become  more  difficult  since 
the  polarization  planes  of  the  694  nm  and  347  nm  laser 
beams  ar'*  normal  to  each  other. 

Finally,  the  PMT  box  can  be  aligned  with  respect  to  the 
optical  axis  of  the  telescope  by  a precision  X-Y-adjust. 
The  coil  (31),  which  is  attached  to  the  photomultiplier 
tube, effects  the  gain  of  the  tube  as  mentioned  later. 

The  operating  voltage  for  the  PMT  is  supplied  by  a highly 
stabilized  Ortec  high  voltagf  power  supply  (32).  The  out- 
put signal  is  fed  to  the  data  acquisition  electronics  via 
a 30  Ohms  terminated  cable. 


Due  to  its  decrease  with  the  squsire  of  distance,  the  return 
signal  traverses  a very  wide  amplitude  range  which  the 
photomultiplier  tube  and  data  acquisition  have  to  cope 
with.  In  extreme  situations,  especially  when  the  scat- 
tering function  itself  grossly  decreases  with  altitude, 
the  PMT  may  be  overloaded  by  the  strong  nearby  signal  and, 
nevertheless,  the  faint  returns  from  greater  distances 
hardly  to  be  measured.  For  this  reason  a device  had  been 
provided  with  the  initial  version  of  the  lidar  system 
which  permits  controlling  the  gain  of  the  PMT  propor- 
tionally to  the  sqiiare  of  the  time  elapsed  from  the  re- 
leasing of  the  Q-switch  pulse.  This  is  done  by  applying  a 
variable  magnetic  field  via  the  coil  (31)  to  the  PMT. 

By  this  field  the  secondary  electrons,  are  more  or  less 
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strongly  deflected  from  the  dynodes  and  so  the  gain  is  re- 
duced. The  current  thi'ough  the  coil  is  generated  in  the 
control  unit  (53)  and  varies  by  a definite  time  function 
which  was  empirically  established  so  as  to  result  in  a 
photomultiplier  gain  increasing  with  the  square  of  time. 
Thur;  the  distance  dependence  of  the  return  signal  is 
eliminated. 

Since  this  magnetic  FMT  gain  control  is  well  operating 
only  over  a gain  range  of  about  2.5  powers  of  ten,  the 
time  range  of  the  device  is  also  limited,  in  view  of  our 
intention  tc  measure  at  higher  levels  of  altitude,  the 
time  raiige  was  fixed  to  range  from  15  to  226  microseconds, 
corresponding  to  an  altitude  range  from  2 to  ^0  km. 

Triggered  by  a reference  pulse  generated  in  time  with  the 
flash-lamp  ignition,  about  1 ms  prior  to  the  v-switch 
pulse,  the  coil  current  goes  to  its  maximum  value,  or  the 
gain  to  its  minimum,  and  remains  then  constant  until 
13  liSfcc  after  releasing  of  the  Q-switch  pulse.  Then  the 
gain  increases  with  the  square  of  time  and  remains  con- 
stant again,  226  psec  after  Q-switching,  at  its  maximum 
value . 

Controlling  amplification  directly  on  the  photomultiplier 
tube  as  opposed  to  a succeeding  electronic  system,  in 
addition  to  simplifying  evaluation  of  the  analog  signal, 
offers  the  advantage  that  an  overloading  of  the  PMT  by 
the  sometimes  very  intense  close-range  signal  is  prevented. 
Of  course,  the  constant  skylight  background  is  also  in- 
creased as  to  interfere  with  the  measuring  signal,  we 

2 

can  work  with  the  t -control  only  at  night. 

4.  Calibration  of  Transmitter  and  Receiver 


In  order  to  be  able  to  compare  lidar  measurements  among 
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each  other  and  also  against  other  date  It  is  necessary 
to  normalize  the  unit.  This  applies  especially  to  two- 
frequency  operation.  The  energy  of  each  individual  laser 
pulse,  both  694  and  347  nm,  is  to  be  measured,  therefore, 
or  monitored.  For  this  purpose,  a small  amount  of  the 
laser  pulre,  reflected  back  from  the  output  window  (20), 
is  fed  to  two  photodiodes  (34)  via  a quartz  fiber  optics 
system  (35).  By  suitable  filters,  it  is  ensured  that  one 
of  the  diodes  is  sensitive  only  to  red  (694  nm)  light, 
the  other  one  to  ultraviolet  (347  nm)  light.  By  this 
reason  only  one  common  output  line  is  necessary  (36). 

A 50  megahertz  A-D  converter  gives  a 4 bit  digital  reading 
of  the  pulse  amplitude. 

Both  diodes  are  calibrated  by  simultaneous  calorimeter 
measurements.  Besides  a conventional  calorimeter  with 
sheet  copper  measuring  and  reference  cones,  a new  volume 
absorbing  calorimeter  is  now  available.  The  laser  energy 
is  absorbed  in  a special  color  glass  plate,  which  avoids 
the  uncontrollable  reflection  losses  of  the  copper  cones. 
Also  recal^bration  is  much  easier  due  to  a special  cali- 
b 'at ion  coil.  The  calorimeter  is  calibrated  by  letting  a 
well-defined  amount  of  electric  energy  flow  for  a short 
time  through  this  coil. 

Since  it  is  necessary  to  measure  the  laser  pulse  snergj* 
at  a certain  minimum  distance  from  the  output  window  (20) 
to  avoid  additional  reflection  back  to  the  photodiodes,  a 
preadjusted  mount  has  been  constructed  to  for  fastening 
the  calorimeter  in  the  proper  position  above  the  laser 
beam  exit  window  (37). 

Figs  III. 3 and  III. 4 show  the  calibration  curves  for  both 
photodiodes.  The  curves  are  linear  but  show  a consider- 
able zero-point  offset.  This  actually  increases  the  accu- 
racy of  the  digital  photodiode  readings  beyond  the  4-bit 
resolution  limit. 

Absolute  calibration  of  the  receiver  is  very  difficult. 
Usually,  lidar  receivers  are  calibrated  relatively  with 


respect  to  the  transmitter  pulse  power  using  retvxrn 
signals  from  an  atmosphere  with  known  backscatter 
function.  This  procedure  is  described  In  detail  in  section 
V of  this  report. 

For  checking  the  linearity  or  possible  gain  changes  of 

the  photomultiplier  tube,  a precision  microsecond  flash- 

lamp  is  available.  This  flashlamp  may  be  triggered  with 

variable  delay,  and  so  it  is  suitable  also  for  checking 

2 

the  performance  of  the  t -control  unit. 

5.  Data  Acquisition 


The  return  signals  may  now  be  recorded  in  three  different 
ways: 

a)  analog  recording  by  means  of  oscilloscope  and 
camera, 

b)  digital  recording  by  means  of  transient  recorder, 
and 

c)  single  photon  counting  with  range  gating. 

5.1.  Analog  Recording 


During  the  first  year  of  operation  of  the  lidar  system, 
only  analog  data  recording  was  possible.  The  signals  de- 
livered by  the  PMT  are  fed  to  a dual  beam  Tektronix  556 
(38)  oscilloscope.  The  traces  on  the  CRT  are  photographed 
either  on  4 x 5"  Polaroid  film  or  70  mm  perforated  nega- 
tive film.  Evaluation  has  to  be  done  manually  in  this  case 
Fig. III. 5 shows  a Polaroid  print  with  examples  of  return 
signals  photographed  directly  from  the  oscilloscope  screen 

5.2.  Digital  Recording  with  Punch  Tape.  Data  Ovitput 


For  digitization  and  storage  of  the  data  a tremsient  re- 
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corder  type  Biomation  8100  (39)  is  provided.  This  recor' 
der  contains  a very  rapid  (100  MHz)  A-D  converter  with 
8 bit  (=.4?6)  amplitude  resolution.  The  memory  capacity 
is  2048  words  d 8 bit.  The  minimum  sample  interval  is^lO 
nanoseconds,  which  results  in  an  altitude  resolution  of 
1.5  m.  The  transient  recorder  is  triggered  in  a manner 
very  similar  to  an  oscilloscope.  In  our  case  it  is  trig- 
gered externally  by  the  trigger  sensor  (8).  Dual  time 
base  operation  allows  recording  of  the  "upper”  part  of 
the  return  signal  with  a slower  sampling  rate  and  thus  ob- 
taining the  exact  amplitude  zero  point,  which  depends  on 
the  varying  sky  background. 

The  stored  data  may  be  output  in  three  different  ways. 

On  the  one  hand,  a reconstructed  analog  signal  output  is 
provided  in  two  forms.  The  Y display  output  presents  this 
analog  signal  repetitively  to  permit  the  signal  to  be  con- 
tinuously displayed  on  a CRT  display  or  oscilloscope.  For 
this  purpose,  a suitable  X output  is  also  provided.  The 
reconstructed  analog  signal  is  also  available  for  output 
to  a strip-chart  recorder  or  Y-T  plotter  (we  are  using  a 
Brush- type  rapid  strip-chart  recorder  (40)  for  this  pur- 
pose, see  Figs  III. 8 - III. 11).  The  plot  time  is  20  sec- 
onds. 

Finally,  via  a digital  output  the  stored  digital  data  may 
be  read,  at  reduced  speed,  directly  into  a computer  in  on- 
line operation,  or  may  be  transferred  via  interface  on 
magnetic  or  punch  tape,  or  the  like.  In  our  case,  data 
output  is  done  in  binary  form  on  8-channel  punch  tape  via 
an  interface  (41)  and  a Facit  tape  punch  (42).  For 
further  processing  the  punch  tape  is  read  into  our  com- 
puter where  the  data  are  normalized  with  respect  to 
transmitter  output  power  and  receiver  sensitivity  (e.g. 

PMT  high  voltage  setting),  and  range  corrected  by  mul- 
tiplying with  the  square  of  distance.  Then  the  converted 
(and  smoothed,  if  necessary)  signal  is  plotted  on  a high- 
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.“:peed  incremental  plotter. 
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Due  to  the  R -law  the  return  signal  covers  a wide  range 
of  amplitude.  As  a consequence  only  signals  from  an  alti- 
tude range  not  exceeding  2 km  can  be  processed  by  the 
transient  recorder  without  a considerable  loss  in  resolu- 
tion. For  this  reason  a gain-switching  amplifier  (43)  is 
available,  which  actually  consists  of  two  amplifiers  with 
gains  of  2.65  and  29.3,  respectively,  operating  in  paral- 
lel. A special  electronic  circuitry  switches  very  rapid- 
ly from  one  amplifier  output  to  the  other.  During  steady 
state  conditions,  the  low  gain  amplifier  is  running. 

When  a laser  pulse  is  initiated,  a trigger  signal  from 
one  of  the  photodiodes  is  received  and  a delay  time  gen- 
erator started.  At  the  end  of  the  presettable  time  delay, 
the  output  is  switched  to  the  high  gain  amplifier,  and 
then  back  again  after  5 milliseconds.  The  switching  is 
very  rapid,  so  that  only  a few  data  words  in  the  tran- 
sient recorder  are  lost. 

5.3.  Examples  of  Data  Records 

Figs  III. 5 - III. 12  show  examples  for  the  different  types 
of  data  recording  and  plotting. 

Fig. III. 5 is  a Polaroid  photograph  of  original  lidar  re- 
turns as  displayed  on  the  oscilloscope  screen.  The  time 
or  altitude  axis  points  up.  The  rotatable,  indexed  sliding 
back  of  the  camera  allows  recording  multiple  images  on 
one  photo.  In  Fig. III. 5 five  backscatter  signals  are  re- 
corded wxth  different  oscilloscope  gain  and  sweep  settings. 
Beginning  with  the  left,  they  are  the  following: 

Signal  No.1:  5 mV/di v. , 10  psec/div.  = 1500  meters/div. , 

" " 2:  5 mV/div.,  5 psec/div.  = 750  meters/div., 

" " 3:  5 mV/div. , 2 psec/div.  ^ 300  meters/div., 

•'  " 4:  20  mV/div.,  2 psec/div.  = 300  meters/div., 

" " 5:  100  mV/div. , 2 psec/div.  = 300  meters/div., 
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All  signals  are  694  nm  returns.  An  aerosol  layer  with  an 
upper  boundary  at  2500  m above  station  level  (3250  m 
above  sea  level)  is  clearly  to  be  recognized,  with  some 
fine  structure  below. 

It  is  obvious  that  the  manual  evaluation  of  the  data, 
using  only  the  oscilloscope  photos,  is  a tedious  and  time- 
consuming  task.  The  Polaroid  photo  Fig. III. 6 represents 
three  pairs  of  simultaneous  694  and  347  nm  return  signals 
as  recorded  with  the  Biomation  8100  transient  recorder. 

The  reconstructed  analog  signals  were  displayed  on  the 
CRT  screen  of  the  oscilloscope  and  photographed,  each 
pair  of  signals  with  coincident  coordinate  aoces.  The 
smaller  signals  are  the  347  nm  (ultraviolet)  signals. 

The  T-axes  point  up,  the  scale  is  Z tisec/div.  or  300  m/ 
div..  The  red  (694  nm)  signals  again  show  considerable 
fina  structure  due  to  aerosol  inhomgenities . The  peak  at 
500  m above  station  level  in  the  0930  CBT  signal  may  be 
generated  by  a transient  small  cloud.  In  contrary,  in  the 
ultraviolet  (347  nm)  signal  the  influence  of  aerosol  is 
much  less  pronounced.  This  is  a feature  common  to  all 
ultraviolet  returns  and  is  surely  caused  by  the  much 
hi<^her  molecular  part  of  the  backscatter  function  at  the 
347  nm  v.avelength  (see  sections  IV  and  V of  this  report). 

Instead  of  taking  photos,  the  reconstructed  analog  signal 
generated  by  the  transient  recorder  may  be  recorded  via 
an  additional  output  connector  by  a fast  Y-T  strip  chart 
recorder,  as  mentioned  above.  We  are  using  a Brush  220 
pressurized  ink  recorder.  The  total  analog  signal  is  out- 
put within  20  seconds.  Using  a chart  speed  of  5 mm/sec, 
this  results  in  an  overall  signal  length  on  the  chart  of 
10  cm.  Figs  III. 7 - III. 10  are  copies  of  original  Brush 
strip  charts,  showing  two  pairs  of  double  frequency  return 
signals.  In  these  cases,  the  gain-switching  amplifier  had 
been  used.  The  gain  switching  is  indicated  by  a very 
rapid  increase  in  signal  amplitude.  The  small  peak  just 
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below  the  switching  point  is  an  artefact  generated  by 
the  switching  circuitry  and  seems  to  be  unavoidable.  It 
is  eliminated  during  computer  processing  of  the  diP'ltal 
data.  Again  the  striking  lack  of  fine  structxire  of  the 
3^7  nm  signals  (besides  the  cloud  peak  at  2500  m a.s.l. 
in  Fig. III. 8),  as  compared  to  the  69A  nm  signals,  is  to  be 
mentioned. 

The  •'shoulder"  in  the  694  nm  signal  curve  Fig.  III. 9 at 
2500  m a.s.l.  (denoted  by  an  arrow),  which  indicates  an 
aeroscl  layer,  is  hardly  to  be  recognized  in  the  corre- 
sponding 347  nm  signal.  Fig. III. 10. 

Fig. III. 11  represents  a pair  of  double  frequency  back- 
scatter  profiles.  The  original  return  signals  are  shown 
in  the  left  part  of  Fig. III. 6.  The  digitized  signals  had 
been  output  on  punch  tape  and  fed  to  the  computer.  Then 
they  were  normalized  and  range-corrected  by  multiplying 
with  the  square  of  distance,  and  plotted  by  our  high- 
speed incremental  plotter. 

Due  to  the  limited  amplitude  resolution  of  the  transient 
recorder,  the  influence  of  digitization  shows  up  ^'ery 
clearly  in  the  upper  parts  of  the  plots,  especially  in  the 
347  nm  profile.  As  mentioned  above,  these  signals  had  been 
obtained  before  the  gain-switching  amplifier  was  avail- 
able. This  amplifier  largely  avoids  the  influence  of  reso- 
lution limitation. 

5.4.  Single  Photon  Counting 


As  mentioned  above,  the  amplitude  of  the  lidar  return 
signals  decreases  rapidly  with  distance  or  altitude.  Due 
to  the  limited  quantum  efficiency  of  the  PMT  cathode, 
beyond  a certain  minimum  distance  depending  on  transmitter 
power  and  receiver  sensitivity,  the  signals  are  more  and 


- 16  - 


more  dissolved  into  individual  photon  peaks.  In  this  case, 
the  sensitivity  and  signal- to-noise  ratio  of  the  receiver 
is  greatly  increased  by  the  range -dependent  phcLon  coun- 
ting technique.  In  oiir  case,  a ten-channel  photon  counter 
is  provided  (44).  By  a special  gating  circuitry,  the  ten 
counters  are  "opened"  in  series  one  after  the  other,  for 
a presettable  gating  time.  The  first  channel  is  opened 
after  an  also  presettable  time  delay  after  the  laser 
pulse  initiation.  An  eleventh  channel,  opened  1 milli- 
second after  the  10th,  serves  as  reference  channel  for 
background  measurement. 

Two  modes  of  operation  are  possible  with  the  counter.  In 
the  one  (single  shot)  mode,  the  stored  photon  counts  are 
output  automatically  after  each  laser  shot  on  punch  tape 
in  BCD  code.  In  the  other  mode  ("adder"  mode),  counts 
generated  by  a presettable  number  of  laser  shots  are 
stored  in  two  additional  storage  registers.  In  detail,  the 
process  is  the  following;  with  n denoting  the  preset 
number  of  single  shot  counts  to  be  added,  during  the  first 
n shots,  the  counts  of  channel  1 and  2 are  added  in  the 
two  storage  registers.  Then  the  content  of  the  registers 
is  output  automatically  on  punch  tape,  and  the  registers 
are  reset.  Subsequently,  during  the  next  n shots,  the 
counts  of  channel  '3  and  4 are  added  in  the  same  way  in  the 
two  registers,  given  out  on  tape,  and  so  on  until  the 
channel  9 and  10  counts  are  added  and  output  on  the  punch 
tape.  This  completes  the  whole  cycle. 

In  order  to  get  single  photon  pulses  with  sufficient 
amplitudes  for  photon  counting  (some  ten  millivolts)  at 
the  auiode  of  the  photomultiplier  tube,  the  tube  has  to  be 
operated  with  relativley  high  voltage  (about  2.2  kilo- 
volts in  our  case).  This  would  imply  a heavy  overload 
of  the  tube  by  the  strong  short-range  signal.  For  this 
reason,  a range  gating  facility  is  provided  with  the  pho- 


ton  coionter.  In  this  mode  of  operation,  dynode  7 of  the 
photomultiplier  is  kept  at  a voltage  which  results  in  a 
nearly  zero  overall  gain  of  the  tube.  Only  dixring  the  gate 
interval  when  the  counters  are  "open”,  this  voltage  is 
changed  to  its  normal  value  by  means  of  an  appropriate 
electronic  circuitry. 

The  range  gating  facility  may  be  useful  also  for  analog 
signal  recording,  e.g.  when  signals  from  higher  altitudes 
are  to  be  recorded  with  good  resolution. 

6.  Fluorescence  Suppression 


In  addition  to  the  coherent  laser  light  the  ruby  also 
emits  incoherent  fluorescent  light  for  about  one  milli- 
second. Although  it  is  more  divergent  and  much  less  inten- 
sive than  the  laser  pulse,  this  fluorescence  would  under 
unfavorable  conditions  (strong  short-range  backscatter) 
superpose  the  faint  return  signals  from  high  altitudes. 

For  this  reason  a mechanical  chopper  (45)  is  provided 
which  blocks  the  beam  path  completely  137  psec  after  Q- 
switching  which  corresponds  to  an  altitude  of  21  km.  The 
chopper  consists  of  a cylinder  with  two  opposite  holes 
and  rotates  wi*h  200  r.p.s.  An  appropriate  electronic 
control  (46)  with  optical  sensors  ensures  exact  synchro- 
nization of  chopper  passage  and  Q-switch  triggering. 

7j A First  Example  for  High  Altitude  Lidar  Measurement 


After  final  completion  of  the  photon  counting  equipment, 
the  weather  conditions  had  been  favorable  for  night-time 
high  altitude  Ixdzr  measurements  for  the  first  time  on 
July  28,  1975. 

At  first,  a series  of  694  nm  analog  signals  were  recorded 
with  the  Biomation  transient  recorder  using  the  reuige 
gating.  The  gate  delay,  i.e.  the  time  lag  between  Q- 
switching  and  the  beginning  of  the  gate  interval,  had  been 


set  to  successively  increasing  values  of  60,  80,  100,  ..., 
300  microseconds,  corresponding  to  altitudes  of  9 to  45 
km  above  station  level  (9.75  to  45.75  km  above  sea  level). 
The  length  of  the  gate  interval  itself  was  30  psec  or 
4.5  km.  This  resulted  in  a set  of  overlapping  partial 
signals,  which  are  shown  as  original  strip  chart  plots  in 
Fig. III. 12.  The  surprisingly  high  amplitude  of  the  signals, 
which  by  far  exceeds  the  time  resolution  limit  of  the 
photon  counter,  clearly  shows  up  in  the  diagram.  Even  from 
more  than  70  km  altitude  a signal  with  non-zero  aimplitude 
is  received,  as  shown  in  the  upper  right  co’mer  of  Fig. 

III. 12. 

The  partial  signals,  although  very  noisy,  can  be  combined 
to  a complete  return  signal.  Fig. III. 13.  The  signal  am- 
plitude decreases  steadily  with  altitude,  but  as  a range- 
correction  clearly  shows,  too  slowly.  The  backscatter 
cross  section  should  decrease  with  altitude  approximately 

proportional  to  the  density  of  the  atmosphere  ( the  den- 
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sity  of  the  standard  atmosphere  decreases  from  .4  kg/m 
It  10  km  a.s.l.  to  3.7  x 10“^  kg/m^  at  40  km). 

Obviously  a nearly  constant  backg.-ound  is  superimposed  to 
the  measuring  signal.  This  background  was  not  generated 
by  residual  skylight,  as  test  measurements  demonstrated. 

It  must  be  assumed,  therefore,  that  the  background  is 
caused  by  ruby  fluorescence.  Since  lidar  returns  from  30 
km  altitude  showed  no  significant  difference  with  and 
without  using  the  fluorescence  suppression,  (cf.  Fig. III. 
17),  it  had  not  been  used  during  the  measurements  de- 
scribed now.  The  first  experiences  seem  to  suggest,  how- 
ever, that  it  is  not  recommended  at  all  to  perform  lidar 
measurements  without  fluoresc'^nce  rejection,  even  under 
apparently  favorable  conditions. 

If  the  probable  fluorescence  background,  which  actually 
doesn't  need  to  be  exactly  constant,  is  subtracted  from 
the  signal,  range-correction  (i.e.  multiplication  with 


the  square  of  distance)  results  in  the  backscatter  pro- 
file shown  in  Fig. III.  14.  Beca'>se  if  the  poor  signal 
resolution  due  to  inherent  noise,  its  fine  structure 
cannot  be  considered  to  be  real. 

For  single  photon  counting  the  signal  amplitude  is  to  be 
reduced  considerably.  For  this  purpose  we  shifted  the  gray 
wedge  from  the  IOO96  to  a 3.4%  transmissivity  position. 

Then  again  a series  of  shots  was  fired  with  the  same  set 
of  range  gate  positions  as  described  above.  For  reducing 
the  statistical  error  ten  shots  were  fired  per  range  gate 
position  and  output  on  punch  tape.  It  is  to  be  borne  in 
mind  that  the  total  gate  interval,  in  our  case  30  psec 
or  4.5  km,  is  divided  by  the  ten  counting  channels  into 
ten  parts  of  equal  length.  This  results  in  an  altitude 
resolution  of  450  m. 

A preliminary'  examination  of  the  results  showed  a sig- 
nificant increase  of  the  coimts  from  channel  1 to  about 
chcinnel  4 or  5,  independent  of  the  gate  delay  setting. 

A somewhat  lower  count  rate  of  the  first  channel  might 
be  due  to  a non-zero  relaxation  time  for  the  gain  of  the 
photomultiplier  tube  to  reach  its  final  value  after  the 
beginning  of  the  gate  interval.  A careful  examination  of 
the  rise  cf  the  PMT  current  after  "opening"  of  the  range 
gat- , when  the  tube  is  exposed  to  a constant  light  source 
Ce.g.  to  the  sky  background),  however,  leads  to  the  con- 
clusion that  thj.s  relaxation  time  doesn't  exceed  a few 
hundred  nanoseconds.  This  results  in  a coiont  reduction 
of  the  first  channel  by  only  a few  percent.  All  other 
channels  should  show  no  differences  in  covint  rate  at  all, 
when  a constant  photon  flux  is  fed  to  the  PMT. 

Maybe  the  mentioned  differences  are  caused  by  improper 
alignment  of  the  threshold  of  the  photon  counter  input 
circuit.  Otherwise,  the  channels  have  to  be  "calibrated" 
by  counting  photon  pulses  generated  by  a constant  light 


source  (sky  background). 

In  Fig. III.  15  the  coionts  of  channels  2 to  10  have  been 
added  and  thus  the  differences  between  the  channels  eli- 
minated. Each  series  of  counts  is  represented  as  2ui  "I”- 
shaped  bar.  Its  vertical  extent  denotes  the  length  of 
the  gate  interval,  the  horizontal  extent  the  statistical 
error,  which  is  equal  to  the  square  root  of  the  sum  of 
the  counts.  The  total  counts  per  gate  interval  steadily 
decrease  with  altitude,  but  tend  to  approach  a non-zero 
background  in  a similar  manner  as  the  analog  signal. 

Fig. III.  13.  .If  this  backgro\md  is  again  subtracted  and 
range-correction  applied,  the  backscatter  profile  shown 
in  Fig. III. 1 6 is  obtained.  The  meaning  of  the  I's  is 
the  same  as  in  Fig. III. 15.  Additionally,  the  density  p 
of  the  steuidard  atmosphere  is  plotted  in  the  diagram  in 
relative  units.  Obviously  the  backscatter  profile  doesn't 
fit  to  the  density  curve  at  mid-altitudes  from  about  20 
to  40  km.  This  is  surely  not  due  to  aerosol  backscatter, 
but  may  rather  be  caused  by  a decrease  of  the  fluores- 
cence background.  On  the  other  hand,  the  step  in  the  back- 
scatter profile  at  30  km  altitude  is  might  be  possibly 
real. 

These  first  high  altitude  lidar  measurements  showed,  on 
the  one  hand,  urgently  that  the  rejection  of  the  ruby 
fluorescence  is  much  more  important  than  it  seemed  to  be 
previously.  On  the  other  hand,  the  extreme  sensitivity 
of  the  photon  counting  method  is  to  be  emphasized,  es- 
pecially if  the  reduced  receiver  sensitivity  due  to  the 
gray  wedge  is  considered.  By  setting  the  gray  wedge  to  a 
higher  transmissivity  position,  a further  increase  of 
the  measurement  accuracy  is  possible  at  higher  altitudes. 

The  high  sensitivity  of  the  photon  counting  method  is 
very  advantageous  to  347  nm  measurements,  which  are  more 
limited  in  range  due  to  the  lower  transmitter  energy  and 
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higher  molecular  extinction  in  that  wavelength.  For 
technical  reasons,  however,  only  694  nm  high  altitude 
data  are  available  up  to  now. 

The  measurements  described  above  didn't  show  any  aerosol 
backscatter  with  certainty.  A night-time  high  altitude 
measurement  performed  in  January  1975,  before  the  in- 
stallation of  photon  counting  and  fluorescence  sup- 
pression, should  be  mentioned  here,  which  revealed  a very 
distinct  backscatter  peak  at  20  km  altitude,  cf.  Fig. III. 
18.  This  peak  is  certainly  due  to  an  aerosol  layer  gener- 
ated by  the  eruptions  of  the  Fuego  volcano  (Guatemala) 
in  October  1974.  Conspicuous  twilight  phenomena  had  been 
additionally  observed  at  that  time  during  dusk  and  dawn. 

f 

IV.  EMPIRI..\L  COMPARISONS  OF  LIDAR  BACKSCATTER  PROFILES 
WITH  CABLE-CAR  SONDE  DATA 

A few  examples  of  empiric  comparisons  between  lidar  back- 
scatter profiles,  on  the  one  hand,  and  results  of  cable- 
car  soundigs,  especially  vertical  profiles  of  Aitken 
nuclei  number  concentration,  on  the  other,  are  presented 
below,  see  Figs  IV. 1 through  IV. 5. 

I In  detail,  the  following  parameters  are  measured  by  cable- 

I car  sondes ; 

> 

I a)  Temperature  (T), 

t 

j b)  Wet-bulb  or  psychrometer  temperature  (T').  This  is 

j the  temperature  of  a ventilated,  wet-tissue  covered 

I thermistor.  It  is  the  more  below  air  temperature  T, 

) the  drier  the  air,  thus  providing  measiire  cf  rela- 

! tive  humidity.  Both  temperature  scales  rise  to  the 

! left. 

i c)  Polar  electric  conductivities  As  mentioned 


I 
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above,  as  a result  of  aerosol  particles  capturing 
small  ions,  conductivity  will  decrease  with  in- 
creasing particle  niunber  concentration.  Thus,  the 
purer  the  air,  the  greater  the  value  of  conductiv- 
ity. Consideration  must  be  given  to  the  fact  that 
high  particle  concentrations  are  measiired  with  a 
lesser  accuracy  than  low  ones.  Conductivity  scales 
also  rise  to  the  left. 

d)  Aitken  nuclei  number  concentration,  N.  A small 

Aitken  nuclei  counter  provides  total  number  of  all 
particles  per  cc,  above  a limiting  value  of  about 
.01  pm  diameter  which  is  determined  by  the  expan- 
sion ratio,  but  otherwise  irrespective  of  size 
distribution.  Again,  the  scale  rises  to  the  left. 

These  profiles  are  compared  against  our  lidar  data  which 
are  represented  by  the  variable  B = PZ  /E  (P  = inc ’ dent 
radiation  power;  Z = distance  or  altitude;  E = transmitter 
pulse  energy)  in  relative  units,  standardized,  however, 
for  each  wavelength  so  as  to  be  comparable  among  each 
other.  As  the  lidar  equation  shows,  term  B or,  respec- 
tively Bp^  for  694  nm  wavelength,  and  Byy  for  347  nm,  is 
proportional  to  the  backscattering  function,  excepting 
extinction  losses.  Because  of  the  short  measuring  distance 
in  this  case  extinction  is  not  assumed  to  have  any  con- 
siderable effect,  at  least  in  the  red.  Unlike  all  other 
parameters  the  B-scale  in  the  graphs  rises  to  the  right. 

In  comparing  B to  N it  must  be  considered  that  the  back- 
scattering  function,  unlike  nuclei  concentration,  is  a 
fvmction  of  particle  size  distribution,  roughly  speaking 
in  the  form  of  preference  to  the  large  particles.  There- 
fore, exact  parallelism  of  the  two  parameters  is  not  to 
be  expected. 

The  profiles  shown  in  Figs  IV. 1-5  were  recorded  11  June 
75  during  a weather  situation  which  is  very  typical  of 
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a moimtaln  valley  opening  into  a plain:  Due  to  heating  of 
the  mountain  slopes  by  insolation  in  the  hours  of  the 
forenoon  air  will  ascend  there.  The  resultant  suction  will 
cause  air  from  the  fore-plain  which  is  more  polluted  and 
frequently  also  more  humid,  to  be  carried  up  the  valley 
more  and  more.  This  process  will  stop  toward  the  evening 
and  be  reversed  in  the  course  of  the  night.  The  flown- 
in  air  is  replaced  by  air  of  a low  aerosol  content, 
flowing  down  from  the  slopes. 

The  Fig. IV. 1 diagram,  at  1020  GET,  still  shows  an  inver- 
sion at  2100  m a.s.l.,  well  defined  by  temperature  in- 
crease and  hiimidity  decrease,  which  had  developed  as  a 
result  of  nocturnal  radiation.  Underneath  we  find,  demon- 
strated by  increased  lidar  backseat ter  and  CN  concen- 
tration plus  decreased  conductivity,  a largely  honge- 
neous  well  mixed  layer  of  an  aerosol  content  higher  than 
that  found  above  the  inversion.  The  stop  in  lidar  back- 
scatter  intensity  is  more  distinct  in  the  red  (Bj^)  than 
in  the  ultraviolet  (B^y). 

This  situation  observed  quite  frequently  may  be  explained 
by  the  considerably  higher  molecular  component  of  the 
backscattering  function  in  the  short  wave  range. 

Closely  below  the  inversion  a maximum  is  found  in  the 
Bj^-curve,  which  is  absent  with  the  Aitken  nuclei.  How- 
ever, it  turns  up  again  as  minimum  in  the  conductivity 
profile.  This  is  assumed  to  be  the  phenomenon  of  pre- 
condensation. Within  the  range  of  the  highest  relative 
humidity  marked  by  the  shortest  distance  between  T and 
T ' , the  particles  have  absorbed  water  to  such  an  extent 
that  increased  backscatter  and  capture  of  small  ions  have 
resulted.  This  interpretation  is  corroborated  by  the 
Fig. IV. 2 diagram,  which  shows  the  conditions  20  minutes 
later.  At  the  inversion  single  cumuli  have  already  devel- 
oped, one  of  which  was  temporarily  located  at  the  zenith 


within  the  field  of  view  of  the  receiver  at  time  of 
the  694  nm  lidar  meastirement  causirig  a correspondingly 
intensive  signal.  Conductivities  here  show  a very  distinct 
minimum,  while  again  no  increased  nuclei  concentration  is 
observed.  With  certainty  the  minimum  of  conductivity  has 
been  caused  by  the  relatively  fast  capture  of  small  ions 
on  cloud  droplets.  Speed  of  capture  of  Aitken  nuclei  on 
cloud  droplets  is  considerably  slower,  hence  a similar 
effect  is  not  observed. 

Some  time  later,  at  1413,  the  inversion  has  disintegrated 
between  2.1  and  2.2  km  a.s.l.,  due  to  increasing  turbu- 
lence. No  structure  whatsoever  is  found  in  lieu  of  it  in 
any  of  the  profiles.  However,  by  the  behavior  of  CN, 

and  air  of  high  aerosol  content  having  flown 
in  from  the  foreplain  can  be  detected  in  the  lower  layers. 
At  its  upper  boundary,  at  1450  m a.s.l.,  development  of  a 
new  inversion  has  commenced.  Different  fine  structures  of 
CN  and  B^^  may  be  caused  by  local  conditions,  inasmuch  as 
we  are  here  dealing  with  a dynamic  process  occurring  in 
the  horizontal  (the  distance  between  lidar  and  cable  car 
sites  is  approximately  10  km). 

About  an  hour  later  (1503  GET,  Fig. IV. 4)  readily  indicated 
by  the  behavior  of  lidar  and  cable-car  data,  the  upper 
boiindary  of  the  air  masses  flowing  up  the  valley  has  as- 
cended to  1550  m a.s.l.  and,  like  the  temperature  in- 
version at  the  same  altitude,  has  become  more  distinct. 
Again  half  an  hour  later  the  upper  aerosol  boundary  has 
arrived  at  1750  m a.s.l.  Below,  the  fine  structure  of 
and  CN  has  weakened.  An  interesting  observation  is  the 
now  existing  minimum  of  particle  concentration  in  the 
zone  of  inversion,  becoming  apparent  in  the  CN  and  con- 
ductivity profiles,  and  Judging  from  the  absence  of  a 
corresponding  minimum  in  the  Bj^  profile  is  apparently 
limited  to  the  small  particles. 


Workers,  such  as  COLLIS  (1970),  KATERLE  (1970)  and 
McCORMICK  (1971),  have  demonstrated  that  it  is  posssible 
to  obtain  aerosol  density  profile  in  atmosphere  using  a 
two-frequency  Ruby-laser,  measiaring  the  backward  scatter- 
ing received  from  different  heights  in  the  atmospheric 
colmn.  They  assumed  the  aerosol  size  to  follow  a Junge- 
distribution  law: 


n(D)s 


dn 

d(lnD) 


= const.  D' 


and  aerosol  particles  to  have  a real  refractive  index. 
Effect  of  humidity  either  on  the  size  distribution  or  on 
refractive  index  was  not  considered. 

The  latest  results  on  aerosol  distribution  by  JUNGE  and 
JAENICKE  (1971)  as  well  as  the  calculations  of  FOITZIK, 
SCHEITHAUER  and  SPANKUCH  (1969)  shew  that  the  aerosol 
size  distribution  can  be  profitably  represented  through 
an  appropriate  combination  of  log-normal  distributions  of 
the  form: 


n{D)  = 


d(lnO) 


fTiTs 


exp  - 


HANEL  (1968,  1970  and  1972)  has  sho™  that  humidity  has 
effect  on  the  size  and  refractive  index  of  aerosol  parti- 
cles. This  in  turn  affects  the  transmissivity  as  well  as 


n»pr>* 
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the  backscattering  function  for  the  transmitted  and  re- 
ceived lidar  signals. 

The  purpose  of  this  chapter  is  to  describe  a mathematical 
structure  to  be  used  for  analyzing  the  signal  data  from 
a double -wavelength  (0.6943  |i  and  0.3472  p)  lidar  sta- 
tioned at  Garmisch,  with  a purpose  to  obtain  the  density 
as  well  as  the  size  distribution  profile  of  aerosol  in 
the  atmosphere.  The  size  distribution  function  of  aerosol 
at  Garmisch,  as  obtained  from  actual  measurements,  has 
been  assumed  a bimodal  log-normal  function  with  the  two 
maxima  at  = 0.38  p and  = 2 p,  and  both  the  loga- 
rithmic standard  deviations,  s^  = Sp  = s = 0.3  (Figs  V.1a 
b,  c).  The  effect  of  humidity  on  size  and  refractive  in- 
dex has  been  dvily  tadcen  into  accovint  and  is  svimmarized  in 
Table  V.1,  and  graphically  represented  in  Fig.V.2,  basing 
our  calculation  on  the  work  of  H^’NEL  (1971). 

Table  V.1 


Relative 

D 

Complex 

Humidity 

°o 

Refr. 

Index 

0 

1 

1.62  - 

0.0200 

i 

0.5 

1.023 

1.605  - 

0.0187 

i 

0.7 

1.070 

1.570  - 

0.0164 

i 

0.8 

1.280 

1.470  - 

0.0097 

i 

0.85 

1.400 

1.420  - 

0.0075 

i 

0.9 

1.520 

1.415  - 

0.0058 

i 

0.92 

1.585 

1.405  - 

0.0052 

i 

0.94 

1.650 

1.400  - 

0.0045 

i 

0.96 

1.84 

1.380  - 

0.0032 

i 

0.98 

2.300 

1.350  - 

0.0015 

i 

1 .00 

1.330  - 

0.0000 

i 

2. 


Basic  Lidar  Equations 


i’he  basic  lidar  equation  for  a backscatter  signal,  con- 
sidering the  scattering  volvime  to  be  a point  soiirce  lo- 
cated at  altitude  H,  expressed  as  power  incident  on  a 
coaxial  receiver,  is  given  by 

W(X.H)=  \ c Q(X)  I F(180®,X.H)  (3) 

^ H 


After  crossover  of  the  laser  pulse  and  receiver  field  of 
view  the  voltaige  at  the  oscilloscope  follows  from  eq.  (3) 
and  is  given  by 

V(X,H)  « Q(X)  A(X)  — F(X,H)  (4) 

H 

with 

AlX)  = y cl  RT{X)  SIX) 


If  the  oscilloscope  voltage  is  normalized  for  the  lidar 

2 

output  energy  Q(\)  and  the  square  of  height  H , expressing 
it  through  a quantity 


V(X.H) 

Q{X) 


Y(X,H) 


then  we  can  write 


H^V(X.H) 

Q(X) 


= A{X)  t2(X,AH)  FIX,H) 


Y(X,H) 


(5) 


The  qxoantities  I and  F are  implicit  functions  of  molecu- 
lar and  aerosol  densities.  Since  we  are  considering  the 
aerosol  distribution  as  double-log-normal,  we  have  two 
niamber  densities  - one  for  each  distribution  - N^(H)  and 
N2(H).  The  total  aerosol  density,  N(H)  = N^(H)  + N2(H). 

3.  Evaluation  of  F(\,H) 


F(X,H)  is  the  total  scattering  function,  the  sum  of  the 
molecular  - and  the  aerosol-scattering  functions: 


F(X,H)  = + fA<X.H) 


(6) 


For  backscattering,  fj^(\,H)  is  given  by 


fM(X.H)- 


9 

4X^ 


_J 

Nm(H) 


LmM{H)2*2. 


3(2^-A) 
6 -7  A 


(7*) 


vnere  A = 0.0350,  for  djry  air,  is  the  depolarization 
factor.  With  nij^(H)  1 , equation  (7*)  approximates  to 


I . 1.0608  (h>m(H)-1)^ 

-N„(H)  ‘ ' 

From  Rawinsonde  values  for  pressure,  temperature  and 
humidity  at  different  altitudes,  the  values  of  mj,j(H) 
(GOODY  1964)  Njjj(H)  and  thereby  fj^jj(X,H)  can  be  calculated. 


i 


V 

i 


I 


The  aerosol  backscattering  function  f.(X,H)  for  a log- 
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normal  distribution  function  dn(D)  given  by  equation  (2) 
is  given  by 


yX,H)  = ydn(D),  J i{a,m^(H)l  dn(D) 


Dso 


with 


4^2  ‘g 


• 

Z^'g' 


> (8) 


C8a) 


'g 


oo 

1 r i(a,mA) 

r'"^i 

^sj  a 

1 

(A 

1 I 

a=o 


da 


(8b) 


various ttjjj,  s and  (corresponding  to 
hiimidity,  g)  have  been  calculated  and  tabulated.  Then, 
for  a double-log-normal  distribution 

F(X,H)  = f„(X,H)  +0  (a„„  s.m  (K))  ■ N,  (H) 

(9') 

* *G2  (“m2-s.mA<H))  ■ Nj(H) 

Or,  in  an  abbreviated  form 


F(X,H)  = f„(X,H)  + 0g,(X,H)  . N,(H)  *03j(X,H)  ■ Nj(H) 


(9) 
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where  , and  <t>Q2  be  calculated  from  the  tabulated 
values  of  i according  to  equation  (8a),  and  fjjj(X,H)  from 
(7). 

4.  Evaluation  of  X^(X,AH) 


for  a given  thickness  AH  of  atmospheric  layer, 
within  which  it  is  assumed  that  the  molecular  and  aerosol 
number  densities  are  constant,  is  given  by 


x^lX.AH)  = exp  - ^2(o,^(H)  ♦ o^(H) ) AhJ 


(10) 


with 


8ti  . ,..i_  8.8869 

T ' 


(11) 


and  0^(H)  for  a log-normal  distribution  given  by 


o^(H) 


j J D^Kg(a,mJ  dn{D) 

D=o 


(12) 


i 

! 


1 
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d,  . ^ K, 
0 G 


(12a) 


and  K 


j Ot  2 

rJaKE(o,mJ  exp-j^;^l  da 
a=o 


(12b) 


various  s and  (corresponding  to 
humidity,  g)  have  been  calciilated  and  tabulated.  Then, 
for  a double-log-nonaal  distribution 


Oa(H)=  dg,  • N,(H).dgj  ■ Nj(H)  03) 

and  T^(\,hH)  will  be  given  by  equation  (10),  with  equa- 
tions (11)  - (13),  represented  as 

x^(X,AH)  = exp  [-2  [o„(X,HI  ♦ dg,(X,H)  • N,(H) 

(14) 

♦dGjIX.H)  . NjIH)]  ah) 


5.  Scheme  of  Calculations 


Equation  (5) 

Y(X,H)  = A(X)  t^CX.AH)  . F(X,H) 


i 
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is  the  basis  of  all  further  calculations.  We  devide  the 
atmospheric  coliunn  above  lidar  into  layers  of  convenient 


I £q.9 

l|Ni.l,I.Ni.2.2  [ 


Y(X.Hi) 


tqS 


Fie.V.,3 


I Eq  9 
. '||NM.Ni2 


Eq.5 

F(y)j;;;^A(X)l 
/ Eq  .14 


Y(X,Hi) 


Fig. V.4 


I ] input 

r ) intermediate 
calculation 


output 


resolution  thickness,  AH,  and  the  lidar  signals  Y(\,H) 
for  two  wavelengths  for  each  height  are  known,  as  also 
fjjj(\,H)  calculated  from  Rawinsonde  data  (7),  then  if  for 
any  height  H.  , the  two  aerosol  densities  N-  . . , N.  „ 
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for  the  bottom  of  the  resolution  thickness  AH,  and  N. . , 

• \ J.  I 

Ni2  for  its  top  known,  X'^(\,AH)  from  equ.  (14)  and 
F(\,H)  from  equ.  (9)  for  both  wavelengths  can  be  calcu- 
lated. V/ith  this  knowledge  and  equ.  (5)  the  system 
constaints  A(\)  for  both  wavelengths  can  be  calculated 
(Fig.V.3). 


Conversely,  if  A(\)‘s  and  2 known  then 

on  calculating  I^(X,AH)  from  equation  (14)  use  of  equation 
(5)  will  give  the  following  two  equations  for  F^(\,H)  in 
terms  of  and  K^2  (Fig.V.4): 


Nn  Nj2  = 


A’^-n-(xJ)^(AH) 


-fjJ(Hi)  (15a) 


krO 


N|,  N|2  = 


YB 


(15b) 


N^.j  and  N^^2  thus  uniquely  determined  from  these  linear 
simultaneous  equations. 

and  N^2  used,  as  described  in  the  above  para- 
graph, to  calculate  ^ and  2'  this  way,  all 

's  and  N^^^'s  - the  aerosol  number  densities  for  fine 
and  coarse  aerosols  distributed  in  a bimodal  log-normal 
distribution  - can  be  obtained  for  (i  = 0...n),  for  a 
value  of  n,  corresponding  to  a height  from  where  the  lidar 
signals  are  still  recvuved  with  a fair  accuracy. 


Another  way  of  determining  the  lidar  system  constants  con- 
sists in  fitting  the  backscatter  signal  to  a Rayleigh- 
backscatter  signal  reduced  by  a constant  extinction  due 
to  aerosol  scattering. 


In  r more  general  fonn  the  basic  lidar  equation  (5)  reads 


Y(X.H)  = A'tX.H„)  t^(X,H-Ho)T^(X.H)lyX,H)  +y(X,Hll 


with 

A'IX.Hq)  =A(X)x^(X,HJ  . 


If  the  aerosol  densities  and  are  known,  the 

extinction-corrected  lidar  system  constants  A*(\,Hq)  can 
be  computed  from 


A'(X.HJ  - 


Y(X,HJ 

■'^(X.H„)[yX.H,)*y(X,Ho)l 


and  the  calculation  of  the  N^'s  is  performed  as  above 
starting  from  height  in  downward  and  upward  direction. 


In  a special  case,  e.g.  above  the  exchange  layer,  aerosol 
scattering  may  be  negligibly  small  compared  to  molecular 
scattering.  Then  there  is  a layer  with  known  aerosol  den- 
sities 


N,(H)  = N2(H)  = 0 if  . 
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For  our  calculations  we  take  AH  = 1 km  and  take  the  first 

two  N-values  from  the  serological  data  at  Garmisch  (740  m 

a.s.l.)  and  Wank-peak  (1780  m a.s.l.)  to  obtain  the  lidar 

R B 

system  constants  A and  A . To  medce  a check  on  these  val- 
ues we  again  calculate  them  for  the  next  adjacent  layer 

using  the  Wank  eind  Zugspitze  (2963  m a.s.l.  ^ aerological 

R 

data  for  the  N-values.  After  finally  establishing  A and 
A as  mean  of  these  two  calculations  the  calculation  pro- 
gram runs  further  to  calculate  the  aerosol  densities  at 
higher  altitudes  as  described  above.  A systematic  flow 
diagram  is  represented  in  Fig.V.5. 

Since  and  N^2  solved  from  linear  equations  (15)  should 
be  positive,  it  sets  a condition  on  the  absolute  values  of 
A^  and  A®.  If  we  solve  (15)  for  and 
conditions: 

N^1>o,  N^2>°»  "then  we  obtain 
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as  a typical  example,  for 

= 4.0624  X 10”^^  X sterad”^ 

^G2  " 5.3653  X 10”^^  X sterad"^ 

Oq-j  = 1.9469  X 10”^^  x sterad”^ 

4>q2  = 3.0871  X 10“'*^  x sterad*'' 

f ^ = 4 X 10”^  X sterad"^ 

= 6 X 10“^  m””*  X sterad””* 

= 4 X 10*^ 

A®  = 4 X 10^ 

- 0.75 

n (x8)2  - 0.5 

k=o  '' 


The  inequality  (16)  becomes 


0.2  < 


Y^-12 

150  Y^- 180 


< 1.7 


-•  Vi»i4m  .*  fc  ^ 


V 

•V. 
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Plotted  on  a graph  in  Fig.V.6,  the  shaded  region  gives 

R B 

the  permissible  values  of  Y and  Y . 


VI.  CALCULATION  OF  RAYLEIGH-  AND  MIE- SCATTERING 

COEFFICIENTS 


Molecular  Scattering 


Radiosonde  sotindings  yield  pressure  p^,  temperature 

u u 

and  relative  humidity  RH^  for  each  baroswitch  click  j. 

From  this  basic  set  of  data  the  corresponding  altitudes 

H.  and  water  vapor  pressures  e . are  calculated. 

J J 

According  to  the  Rayleigh  scattering  theory  (GOODY  1964) 
the  molecular  volume  scattering  coefficient  is  given  by 


. 32ti’  ,-4  (mj(X)-l)^  3(2.A) 

° ~ ^ — Nj rw 


p^[mb] 

with  molecular  density  N.  = 0.72435  • 10  ^ -z  ^ 

J T,  [K] 


m 


-3 


The  index  of  refraction  is  related  to  pressvire,  tempera- 
ture amd  water  vapor  pressure  by 


(mj(X)-D  -10®  = A(X)  -6(mj(X)-l) 

with  A (0.6943p)  = 78.39, 

A (0.3472p)  = 81.38, 


(p  . is  measured  in  mb  and  T.  in  K) 
J J 
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and  the  water  vapor  pressure  correction 


6 (m.(\)-1) -10* 


e . [mb] 

12.5  ^ = 0.6943P  , 

T^[K] 


e.[mb] 


The  mean  depolarization  factor  for  air  is  A = 0.035. 


Thus  we  obtain  the  molecular  volume  scattering  coeffi- 
cients 


.(0.6943u)m'''  = 1.2C-  10~^(1-0.16  , 

p t[K  ] 


0„  .(0.3472p)m"''  = 2,21  • 10~^(1-0.1'^  -)^  . 

M.O  p 


The  transmissivity  Ij(4)  is  computed  from  the  scattering 
coefficients  by 


Xj(\)  , exp  (-1  a«J  , . 0. 


Two  tables  are  built  up,  one  containing  altitude,  tempera- 
ture, relative  humidity,  Rayleigh  scattering  coefficients 
and  transmissivities  for  each  baroswitch  click  (table 
RAY1 ) (Tab.VI.1),  another  one  containing  H. , RH- , 0^  . 
and  0^  i ^ series  of 

fixed  levels  H.  , e.g.  from  1000  m a.s.l.  to  3000  ni  a.s.l. 
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Tab.  VI.  1 


Rayleigh  scattering  coefficients  and  transmissivities 


26.02.75/09.36 


n 

T 

RF 

sicha-rot 

TAU-ROT 

SIGHA-BLAU 

TAU-BLAU 

U0735 

.'4.2 

085.0 

*65  36E-05 

.77668E-06 

0Q780 

-ns.- 

072.0 

•66866F.05 

.99976 

.77370E-  4 

.99581 

U0882 

-06.7 

076.3 

*662s6E-05 

.99935 

.76371E-06 

.98880 

P0971 

.05.0 

077.4 

•63835E.05 

.99896 

.75594E-P6 

.98215 

01061 

.05.0 

075.9 

•63334E.05 

.99857 

.74730E-06 

.97556 

01153 

.05.1 

072.4 

.62859E.05 

.99817 

.73911E.06 

.968<'7 

01263 

.05.': 

068.0 

•62260E.05 

.99771 

.72642E.04 

.96123 

01352 

.06.7 

064.6 

•61716E.05 

.99734 

.71938E.06 

.95507 

01A52 

.05.0 

062.6 

.61239E-05 

.99693 

.71117E-06 

*96836 

01543 

.05.1 

059.0 

•60788E-05 

.99656 

.70337E-06 

.96227 

01636 

-05.1 

053.6 

.60313E.05 

.99619 

.69518E.06 

.93623 

01753 

-05.0 

047.3 

•39696E.05 

.99572 

.68454E-06 

.92875 

01843 

-03.6 

041.6 

.39063E-05 

.99537 

.67326E-P4 

.92311 

01940 

-03.6 

038.2 

•38569E-05 

,99500 

.6651PE-06 

.91720 

02033 

• 06.4 

038.1 

«38228E-05 

,99465 

.65922E-P4 

*91161 

02131 

-05.1 

037.5 

.37861E-05 

.99428 

.65287E-06 

.90576 

02251 

-05.9 

036.8 

.37396E.05 

.99383 

.64482E.04 

.89881 

02346 

.06.8 

035.8 

.37065E-05 

.99348 

.63915E-04 

.89333 

02443 

-07.5 

035.6 

.36711E.05 

.99313 

.63305E-O4 

*88780 

02535 

-08.4 

035.9 

•36399E-05 

.99279 

.62766E-P6 

.68275 

02634 

-09 . 0 

038.0 

«36019E.05 

.99244 

.6211PE-P4 

.87730 

02733 

-09.9 

040.5 

.35689E-05 

.99209 

.61542E.06 

.87203 

02828 

-10.4 

045.6 

.35314E-05 

.99176 

.6P894E-06 

.86700 

02929 

.10.9 

045.7 

.34908E-05 

.99141 

.6P195E-04 

.86172 

03026 

-11.8 

044.3 

.34595E-05 

.99107 

.59655E-04 

.85676 

03125 

-11.9 

043.0 

.34161E-05 

.99074 

. 5 8906E  —06 

.85180 

03246 

.10.9 

041.3 

•33693E-05 

.99034 

.57754E-06 

.86583 

03342 

-10.7 

041.0 

.33059E-05 

.99002 

.57007E-P6 

.86122 

03-39 

-10.9 

041.7 

.32662E-05 

.98971 

.56322E-06 

.83666 

03543 

.11.6 

042.7 

.32311E-05 

.98930 

.55716E-04 

.83181 

03636 

.12.1 

043.1 

*31989E-05 

.9B90ti 

.55161E-P6 

.62756 

03766 

-12.3 

043.8 

.31466E-05 

.98860 

.54259E-06 

.82173 

03868 

-12.8 

043.6 

.31116E-05 

.98837 

.53656E-06 

.81725 

03971 

-13.5 

044.2 

•30781E-05 

.98805 

.53079E-P4 

.81280 

04075 

-14.0 

045.2 

.3a41lE-05 

,98774 

.5266PE-06 

*80837 

04174 

.14.5 

045.0 

*30qeiE-05 

.98744 

*5lB7lE.06 

.80622 

06281 

-15.4 

044.7 

.29756E-05 

.98713 

.51311E-06 

.79986 

06382 

-15.6 

045.1 

•29390E-05 

.98684 

.5P68PE-06 

.79576 

06485 

.16.6 

046.9 

.29102E-05 

.98654 

.5P182E-P6 

.79165 

04589 

-17.0 

046.5 

.2876BE.05 

.98625 

.69573E-P6 

.78760 

04700 

-17.3 

046.0 

•28354E-05 

.98594 

.68892E-06 

.78331 

04827 

-17.7 

045,3 

•27923E-05 

.98559 

.6815PE-P6 

.77855 

04935 

— 18 . 6 

045.0 

•27625E.05 

.98529 

.67636E-06 

.77656 

05037 

-19.2 

042.9 

.27307E-05 

.98502 

.67086E-P6 

.77083 

05148 

-20.2 

040.4 

«27019E-05 

.98473 

.66590E-P6 

*76688 

05245 

-20.5 

039.6 

•26693E-05 

.98447 

.66028E-P6 

.76346 

05380 

-21.6 

038.8 

.26369E.05 

.98412 

.45434E-06 

.75879 

05487 

-22.3 

037.2 

.26011E-05 

.98385 

.44851E-P6 

.75515 
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in  steps  of  50  m (table  RAY2)  (see  the  first  four  col\amns 
of  Tab.VII.l).  Table  RAY2  is  stored  on  magnetic  tape  or 
disk. 

2.  Aerosol  Scattering 

Originally  the  aerosol  extinction  efficiencies  Kg  (c.f. 
V.4)  were  calculated  after  PENNDORF  and  DEIRMENDJIAN 
formulae  given  in  HiiNEL  (1970).  The  computation  of  the 
aerosol  backscattering  function  f^  requires  the  evaluation 
of  the  MIE-function  for  backscattering 


i(a,m,180°)  = I I (n+1/2)  • (-1)"+^  • (a  -b  ) \ 
n=l  “ 


where  a is  the  size  parameter,  m the  complex  index  of 

refraction  and  a^  and 
n 

coefficients  given  by 


refraction  and  a„  and  are  the  complex  HIE-acattering 


a„(a.m) 


("rn  a) ' g"  ■ gn-i 


(mA„  ♦ f • gn  - gn-1 


The  fxinction  Aj^(ma)  can  be  computed  using  the  recurrence 
relation 


An(ma)  = - =7:  + — ;r 

mci  n . / _\ 

moi 


1 


j 


and 


A^Cma)  = cot(ma). 

The  fvinction  is  also  given  by  recurrence: 


|„(al  - . 


and  = cosa  -i  • sintt  , 

5 q(0L)  = sina  +i  • cosa  • 


As  the  extinction  efficiency  can  also  be  expressed  in 
terms  of  the  ftinctions  a^^  and  b^^  by 


Kg(a>m) 


^ n^l  ‘’n>  ’ 


extinction  efficiency  and  backscattering  MIE-function 
will  be  evaluated  and  tabxilated  in  the  same  nm  for  a 
given  set  of  size  parameters  and  index  of  refraction. 
Numerical  integration  of  these  tabulated  values  over  the 
assumed  log-normal  particle  size  distribution  then  yields 
the  required  collective  values  and  i^  for  given  re- 
fractive index  m and  a log-noiuial  size  distribution  char- 
acterized by  its  mean  size  parameter  ttjjj  and  standard  de- 
viation s (Tab. VI. 2). 
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VII.  THE  EVALUATION  ROUTINE  - COMPUTATIONAL  FLOW  AND 
APPLICATION  TO  ABSOLUTE  LIDAR  CALIBRATION,  AER- 
OSOL DENSITY  AND  SIZE  DISTRIBUTION  PROFILES 

1 . Data  Processing 

Data  processing  occurs  in  two  steps.  In  the  first  step 
the  digital  backscatter  data  edited  on  paper  tape  by  the 
fast  transient  recorder  have  to  be  normalized  and  correc- 
ted for  the  square  of  height,  smoothed  and  plotted  for 
control  and  qualitative  evaluation.  The  backscatter  data 
are  stored  on  magnetic  tape  or  disk  in  two  tables,  one 
for  each  lidar  frequency. 

In  a second  step  lidar  data  and  Rayleigh  scattering  data 
are  combined  into  one  table,  calibrated  and  evaluated 
according  to  the  iteration  method  described  in  section 
V.5  in  order  to  obtain  the  aerosol  profiles. 

a)  Normalization  of  the  backscatter  data 

The  output  of  the  transient  recorder  consists  of  digital 
backscatter  signals  Y'.  (\)  for  2000  sample  intervals  j per 
frequency  K,  divided  into  two  ranges  of  different  amplifi- 
cation and  two  ranges  of  different  time  base. 

The  second  time  base  - for  the  last  200  sample  intervals  - 
f is  comparatively  long  (500ns  = 75m)  and  it  is  assiimed 

' that  the  mean  signal  contained  in  the  last  50  sample  in- 

tervals represents  the  skylight  background  and  determines 
the  zeroline  of  the  second  amplification  range.  The  con- 
dition of  steadiness  at  the  switching  point  betweei  the 
two  amplification  ranges  together  with  the  known  ampli- 
fication ratio  enables  the  computation  of  the  zeroline  in 
the  first  amplification  range. 

t 

'[ 

i 
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The  height-corrected,  normalized  backscatter  signal  is 
given  by 


Yj  « Q-Yj  *(Hf  + b)/E, 


where  ; 


E : 
Q : 


b ; 


background-corrected  signal, 
mean  altitude  corresponding  to  sample 
interval  o, 
emitted  laser  energy, 

normalization  constant  containing  photo- 
multiplier h.v.,  amplification  factor, 
transient  recorder  input  sensitivity, 
defocussing  correction  for  signals  received 
from  lower  altitudes. 


In  order  to  obtain  a reasonable  normalized  backscatter 
signal  and  calibrated  aerosol  profiles,  several  para- 
meters may  have  to  be  varied  during  the  calculations  and 
the  effect  of  these  variations  at  the  different  stages  of 
the  data  processing  shculd  be  visualized  immediately.  To 
this  purpose  a flexible  operation  system  was  developed 
for  our  minicomputer,  using  the  graphical  display  teimi- 
nal,  fast  incremental  plotter  and  fast  line-printer.  The 
following  flowcharts  (Fig.YII.1)  indicate  the  regular  flow 
of  calculations.  Each  operation  in  this  sequence  may  be 
restarted  with  changed  parameters  and  is  controlled 
through  instant  graphical  output. 

b)  Iterative  computation  of  aerosol  density  and  size 

distribution  (cf.V.5). 

The  main  iteration  program  is  operating  on  the  complete 
table  containing  altitude  H. , relative  humidity  RH- , mo- 

1 R ^ B 

lecular  scattering  coefficients  0^^  ^ for  red  and  0^  ^ 


V 


- 4fj  - 


for  UV  lidar  frequencies  and  height-corrected  lidar  back- 
scatter  signals  and  Yj^  for  i = ^•••^max'  P^irameters  of 
the  calculations  are: 

- the  mean  size  parameters  CL^  ^ and  2 stand- 

ard deviations  s^  sind  S2  of  the  two  log-normal  distrib- 
utions characterizing  the  aerosol  size  distribution, 

R B 

the  calibration  constants  A and  A (or  more  general 
A*^'  and  A®'  (cf.V.5)  )• 

For  our  first  calculations  we  choose  size  parameters 
(X  = TID/X  corresponding  to  particle  disimeters  = C.38p 
and  D2  = 2p  and  standard  deviations  s^  = S2  = 0.3,  whicn 
fit  well  to  the  size  distributions  measured  at  our  three 
stations  (see  also  BARTUSEK  et.al.  (1970)  ). 

As  the  index  of  refraction  and  therefore  the  extinction 
efficiency  dG  (cf.V.4)  and  scattering  function  <t>(cl.V.3) 
depend  on  relative  hiimidity , a table  of  dG ' s and  <l> » s as 
functions  of  relative  humidity  has  to  be  provided  for  the 
calculations.  For  the  moment  we  only  use  values  of  dG  and 
calculated  for  one  complex  index  of  refraction  m = 1.62 

- 0.02i,  supposed  to  be  valid  for  a relative  humidity 
rsinge  between  0 - 50%  . 

Finally  the  starting  values  ^ = N^(Hq)  and  ^2  ^ = 

^*2^^o^  have  to  be  fixed  either  by  measurement  at  one  of 
our  mountain  stations  or  by  the  assumption  of  a pure 
Rayleigh  scattering  case  at  H = H^. 

The  flowchart  of  the  main  program  is  given  below  (Fig. VII. 
2).  The  program  calls  several  subroutines.  One  of  them  is 

p 

intended  to  look  up  the  values  for  dG(RH^,(lt^,  s.j ; , 
dG(RH^,a^,  S2),  dG(RH^,a®,  s^ ) and  dG(RH^,  a®,  S2)  and 
0(RH^,a^,  s^),  0(RH.  , a^,  S2),  <t>(RH.  , a®,  s^)  and 
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<t>(RH^,  Ct®.  S2)  in  a precomputed  table  for  a given  set  of 
humidity  intervals.  The  other  subroutines  are  function 
subroutines  to  evaluate  the  aerosol  baokscattering  func- 
tions f^  ^ as  ftinctions  of  Y^,  0^,  and 

^2  i-1  "the  formulae 


Yj 


At?  . x2  . 

i^ti  M,i 


* f, 


where 


3_ 

8Tl 


xf  =xf_,  exp(-20iAHj) 


and 


dG,  + N2,i.idG2. 


The  output  of  the  iteration  program  includes  computer 
plots  of  the  -two  density  profiles  N<j(H)  and  N2(H)  for 
fine  auid  coarse  aerosol  (Fig. VII. 4)  and  a tabulation  of 
the  complete  data  set  (Tab.VII.1),  which  is  the  basis  of 
the  iteration  calculations. 

2.  A First  Example  for  the  Application  of  the  Evalua- 
tion Routine 

The  example  is  taken  from  a series  of  lidar  firings  made 
on  February  26,  1975  at  0940  GET,  It  coincides  with 
a radiosonde  ascent  (Tab.VII.1,  first  part), 
data  transmission  from  the  Zugspitze  cable-car  tele- 
metry system  (Fig. VII. 4), 

simultaneous  measurements  of  the  aerosol  size  dis- 
tribution at  the  three  stations  (Fig.VII.3a,b,c). 
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Tab. VII. 1 


Particle  density  profiles  from  lidar  backscatter  signals 


26.02.75/09.40 


H,  m 

RH,  96 

°m*  ^ 

°m’  “ 

yR 

N.J  ,part/m^ 

Np.part/m 

0O755 

Cb  ,2 

.99979E-05 

.77558E-09 

17.13 

29.85 

0,10E-^P9 

0.16E-f06 

cuo  0 

076.9 

.93675E-05 

.75318E-09 

35.99 

95.i-3 

0,37E-^C9 

0.76E-^06 

1050 

076.  1 

.93398E.05 

.79890E-09 

39.81 

99.11 

0.38E-t^09 

0.6lE-f06 

nil  0 

079.9 

.93139E-05 

.79389E-09 

33.68 

91.51 

0.35E+09 

0.75E+06 

U1150 

072.5 

.92879E-05 

.73936E-09 

3 .98 

38.77 

0.32E-I-09 

0.70E+06 

012  0 

07  .5 

.92599E-05 

.73953E-09 

28. 6>/ 

38.99 

0.34E-«-09 

0.27E•^06 

•1250 

068.5 

.92313E-05 

.72968E-09 

28.52 

36.39 

0.32E^09 

0.52E-^06 

U13  0 

066  « 6 

.92023E-C5 

.72967E-09 

26.96 

32.83 

0.26E-I-09 

0.79E.f06 

•J1350 

069.7 

.91729E-05 

.71962E-09 

26.2' 

32.90 

G.27E-f09 

0 . 68E  4-06 

' v'14  0 

C63.7 

.91987E-05 

.71599E-09 

25.87 

29,65 

0.22E•^09 

O.llE-^07 

■H'iSO 

062.8 

.91298E-05 

.71131E-09 

29.71 

27,59 

O.lBE-^09 

0.12E■^07 

'•'15  0 

■:t>  .8 

.910  IE-O5 

.707  5E-09 

29.08 

26.63 

0.17E+09 

0.12E+07 

91550 

058.6 

.90753E-C5 

.702 /6E-09 

22.50 

2 5. 05 

0.15E-809 

0.12E•^07 

‘•l^  0 

055.7 

.90996E.05 

.69833E-09 

17.02 

23.23 

0.12E+O9 

0.63E-fC6 

‘..11650 

052.8 

.90237E-05 

.6<J3a7E-09 

13.52 

.1.16 

0,86E-*-08 

O.4IE4O6 

'U7  0 

05  .2 

.39979E-O5 

.68933E-09 

12.39 

19.60 

0.55E■^08 

0 . 48E  -406 

01750 

097.5 

.39711E-05 

.68979E-09 

1 •■.28 

19.  18 

0,57E+06 

0. 19E406 

018  0 

099.3 

.39355E-05 

.67865E-09 

09.98 

13.71 

0.55E■^08 

0.86E405 

O1850 

091.9 

.39010L-C5 

.67269E-09 

09.26 

18.39 

0.54E-t-08 

0. 798405 

019  0 

039.6 

.38765E-05 

.66897E-09 

09.03 

17,0', 

0,25E-t-08 

0.27E406 

01950 

038.2 

.38532E-05 

. 66996E  — 09 

08. 8 f- 

16.69 

0.21E-»^0e 

0.29E406 

020  0 

•■'38.1 

.3839tE-05 

.66128E-09 

08. 7U 

16.97 

0.21E+08 

0.27E406 

U2050 

038.0 

.38169E-05 

.65810E-09 

08.23 

17,95 

0.58E+08 
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Although  lidar  data  were  recorded  up  to  a maximum  height 
of  6750  m a.s.l.  and  radiosonde  data  are  available  up  to 
29  km,  the  evaluation  was  limited  to  a maxi mum  height  of 
3000  m a.s.l..  This  is  the  range  where  detailed  compari- 
son with  aerological  and  aerosol  measurements  is  possible. 

This  comparison  is  shown  in  Fig. VII. 4.  The  aerosol  density 
profiles  fit  fairly  well  to  the  particle  concentrations 
derived  from  the  impactor  measuirements  carried  out  at 
Weink  peak  (denoted  by  and  ¥2)  and  at  Zugspitze  (de- 
noted by  and  Z2).  This  is  the  result  of  the  calibra- 
tion routine,  the  aim  of  which  is  to  find  a consistent 

set  of  pairameters  for  both  lidar  frequencies,  viz.  the 

R R 

lidar  system  constants  A and  A , the  zerolines  for  both 
amplification  ranges  and  the  defocussing  correction  con- 
stants. Fitting  the  computed  aerosol  density  profiles  to 
the  measured  values  is  not  yet  carried  out  completely 
automatically,  because  the  change  of  a parameter  for  one 
lidar  frequency  affects  both  aerosol  density  profiles 
(the  backscatter  intensity  of  light  of  any  frequency  is 
due  to  coarse  and  fine  aerosol)  and  may  have  to  be  com- 
pensated pai’tly  by  a oriange  of  a pai'aiiieter  for  the  other 
frequency  etc.  in  order  to  reach  consistency.  The  effect 
of  this  interdependence  of  both  profiles  N.|(H)  and  N2(H) 
Ccin  be  seen  in  Fig. VII. 4.  An  ii:crease  in  the  fine  aerosol 
concentration  is  accompanied  ratner  frequently  by  a de- 
crease in  the  coarse  aerosol  concentration  or  vice  versa. 
Certainly  this  is  not  always  true  in  reality  sind  may  be 
the  result  of  a fluctuation  in  one  of  the  backscatter 
signals . 

On  the  other  hand  a simultaneous  decrease  in  both  fine 
and  coarse  aerosol  densities  indicates  a real  effect  as 
can  be  seen  from  a comparison  between  the  atmospheric 
conductivity  profiles  and  the  aerosol  profiles  in  Fig. 

VII. 4.  In  this  example  the  conductivity  profiles  together 


I 


with  the  temperature  profile  suggest  the  presence  of 
roughly  fovir  layers  between  1000  m and  3000  m a.s.l.. 

The  haze  layer  (1000  - 1500  m)  is  characterized  by  a high 
concentration  of  both  fine  and  coarse  aerosol . Above  the 
haze  layer  there  is  a layer  with  increasing  atmospheric 
conductivity.  This  means  a decrease  of  particle  number 
concentration.  Actually  the  aerosol  concentration  pro- 
files obtained  from  lidar  backscatter  data  show  this  be- 
havior. Above  the  inversion  at  1900  m the  conductivity 
decreases  again.  Possibly  the  increase  if  coarse  aerosol 
concentration  at  this  altitude  together  with  the  fine 
aerosol  density  remaining  constant  may  be  interpreted  as 
a precondensation  effect  within  this  layer.  In  the  fourth 
layer  (2200  - 2800  m)  the  coarse  aerosol  concentration 
decreases  to  a value  of  the  same  order  as  the  one  meas- 
ured at  Zugspitze  peak.  The  absence  of  fine  aerosol  data 
simply  means  that  in  the  course  of  the  evaluation  routine 
negative  particle  concentrations  occurred  probably  due  to 
an  overestimation  of  the  347  nm  backscatter  signal  back- 
ground. At  this  stage  the  iteration  computations  should 
be  restarted  with  changed  parameters.  As  mentioned  above, 
the  then  resulting  Incx'ease  in  the  fine  aerOSOl  COnOcxi— 
tration  would  necessarily  be  accompanied  by  a decrease 
in  the  coarse  aerosol  density. 
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Definition  of  Symbols 

Note:  The  superscript  R or  B above  a symbol  indicates 

that  it  corresponds  to  the  wavelength  0.6943ji 
or  0.3472ti,  respectively. 


Symbol 


Definition 


D 

D, 


m 


g 

H 

Kp 

K. 


1 

m, 


Constant  for  lidar  system 

Extinction  coefficient  of  a single  particle, 

- 

Elxtinction  coefficient  of  log-normal-distri- 
buted particle  collective,  normalized  for  a 
single  particle,  = Kq 

Diameter 

Median  diameter  of  a log-normal  distribution 

Voliime  scattering  function  for  aerosol 

Volume  scattering  function  for  atmospheric 
molecules 

Relative  humidity 

Height  above  lidar  station 

Extinction  cross-section  of  a particle 

Extinction  cross-section  of  a log-normal 
distribution  of  particle  collective,  normalized 
for  a single  particle;  equation  (12b) 

Mie-function  for  backscattering 

Mie-function  of  a log-nomal  distribution  of 
particle  collective,  normalized  for  a single 
particle,  equation  (8b) 

Complex  refractive  index  of  aerosol  and  air, 
respectively 
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Definition  of  Symbols  (cont) 


n(D) 

“h 

N^,  N2 


P 

Q 

r 

R 

S 

s 

T 

W 

a 

\ 


X(AH) 

4>g 


Aerosol  size  distribution  fxmction 

Molecular  niomber  density 

Aerosol  number  densities  for  the  fine  and 
coarse  aerosols,  respectively,  of  a bimodal 
log-normal  distribution 

Pressure 

Lidar  pulse  energy 
Aerosol  radius 

Load  resistance  of  the  photomultiplier 

Spectral  sensitivity  of  the  photomultiplier 

Logarithmic  standard  deviation  of  log-normal 
particle  size  distribution 

System  optical  efficiency 

Lidar  power  signal 

Scattering  parameter  = 

K 

Wavelength 

Molecular  eind  aerosol  extinction  coefficients 
respectively 

Transmissivity  of  atmospheric  layer  thickness 
AH,  equation  (10) 

Volume  backscattering  function  for  a log-nor- 
mal distribution  of  particles,  normalized  for 
a single  particle,  equations  (8a)  and  (8b) 

Effective  area  of  the  receiver 


1 


1 


■y> 


'l/  Vr 
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VIII.  CONCLUSIONS 


1 .  Lidar  System 


The  technical  development  of  the  lidar  system  may  now  be 
considered  as  to  be  essentially  finished.  In  its  present 
layout  the  system  meets  all  requirements  to  be  made.  The 
reliability  of  its  operation,  however,  should  be  further 
improved  by  providing  more  spare  parts,  especially  opti- 
cal components  (ruby,  mirrors,  glass  elates,  frequency 
doubler  etc.)  and  some  electronic  components. 

Future  developments  in  electronics (e.g.  A-D  converters 
with  increased  resolution)  and  optics  (more  efficient 
and  notwithstanding  durable  frequency  doubling  materials) 
should  be  considered  and  utilized. 

2.  Data  Processing 


Our  computer  system,  after  all  programs  being  essentially 
completed,  is  now  fully  capable  of  performing  all  data 
evaluation  work  as  well  as  all  the  extensive  calculations 
described  in  sections  V to  VII  of  this  report.  Neverthe- 
less, additional  on-line,  real-time  data  processing  would 
be  desirable,  especially  in  connection  with  photon  coun- 
ting. In  this  case,  it  is  more  difficult  to  get  a first 
survey  of  the  results  immediately  after  the  measurement. 

3.  Theoretical  Procedures 


The  fundamentals  of  the  theoretical  evaluation  and  cali- 
bration procedures  may  be  considered  as  completed.  With 
the  experiences  gained  at  the  evaluation  of  the  systema- 
tic calibration  measurements  to  be  performed  in  the  near 
future,  some  changes  and  improvements  could  turn  out  to 


be  necessary.  The  main  difficulty  in  the  procedure  is  the 
somewhat  unsatisfactory  resolution  of  the  347  nm  signals, 
even  if,  due  to  minor  defects  of  some  optical  parts,  the 
transmitter  output  power  in  that  wavelength  was  reduced 
at  the  time  the  first  calibration  measurements  were  per- 
formed. With  the  present  development  of  new,  more  effi- 
cient types  of  frequency  doublers,  improvements  in  this 
field  are  to  be  expected  in  the  future. 

4.  Calibration  of  the  Lidar  System 


The  first  results  available  of  absolute  calibration  meas- 
urements using  aerosol  data  from  the  mountain  stations 
indicate  the  procedure  to  be  applicable  in  principle.  It 
has  to  be  improved,  however,  in  some  details,  after  much 
more  empirical  data  for  calibration  have  been  collected. 

IX.  FUTURE  PLANS 

a)  As  the  lidar  system  is  now  essentially  completed, 
only  a few  technical  improvements  are  to  be  carried 
out  in  the  future: 

on-line  data  processing,  new  A-D  converters,  using 
more  efficient  frequency  doublers. 

b)  The  reliability  of  the  system  is  to  be  further 
improved,  among  other  things  by  providing  additional 
spare  parts. 

c)  systematic  calibration  measurements  are  to  be 
carried  out,  together  with  simultaneously  obtaining 
of  aerosol  and  aerological  data  by  the  mountain 
stations,  by  the  cable-car  equipment  and  by  means 
of  the  radiosondes. 

d)  With  the  calibration  of  the  system  completed,  it  is 
to  be  used  for  practical  work  in  remote  sensiiig  of 
aerosols,  including  systematic  night-time  high  al- 
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titude  measurements  with  photon  counting.  In  these 
cases,  simultaneous  night-time  radiosonde  ascents 
are  possible  and  necessary,  since  the  routine 
weather  service  data  are  insufficient  for  this 
pm'pose . 

e)  The  theoretical  procedures  worked  out  in  the  past 
are  to  be  applied  for  evaluation  of  data,  and  to 
be  improved  if  necessary.  The  tables  of  Mie-func- 
tions  should  be  extented  by  calculating  the  func- 
tions for  more  values  of  the  complex  index  of  re- 
fraction. The  influence  of  humidity  is  to  be  con- 
sidered. Finally,  high  altitude  measurements  with 
photon  counting  must  be  included. 

f)  A second,  mobile  lidar  system  is  under  construction 
and  will  be  available  in  the  beginning  of  1976. 

This  system,  besides  a ruby  laser  with  frequency 
doubler,  will  be  equipped  with  an  additional  fre- 
quency doubled  Neodjrmium  glass  laser  allowing 
firings  in  553  nm  wavelength.  With  this  mobile 
three  frequency  lidar,  especially  in  connection 
virith  the  first  system,  vertical  and  horizontal 
variations  of  aerosol  constitution  as  well  as 
processes  of  cloud  formation  and  disintegration 
are  to  be  studied.  The  third  wavelength  will  in- 
volve considerable  modifications  of  the  mathemati- 
cal procedures. 


X. 
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